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ABSTRACT 

We extend a semiclassical model of transfer reactions to the 
case in which one of the collision partners is a vibrational nucleus. 
The model is applied to one- and two-proton stripping reactions in 
the 37 C1 + 98 Mo system, for which a rapid transition from normal 
to anomalous slope in the two proton transfer reaction at energies 
around the Coulomb barrier is experimentally observed. This 
behavior is satisfactorily reproduced by the present extension of 
the model. 

PACS number(s): 25.70.Hi, 24.10.Ht 
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I. INTRODUCTION 



The standard treatment of subbarrier nuclear transfer reactions ffl con- 
siders this process as tunneling of a particle from the potential well created 
by the donor core to the one of the acceptor core. For a given scattering 
angle the tunneling is dominated by the contribution from the associated 
distance of closest approach in the classical Rutherford trajectory, Drum, 
and the transfer probability is 



where fi and B e ff are the reduced mass and the effective barrier height to be 
traversed by the transferred particle, respectively. When, as customary, the 
transfer probability is presented as a function of DRuth in a semilogarithmic 
plot, this model predicts a straight line with an energy independent slope, 
being the slope for two-nucleon transfer approximately twice that for one- 
nucleon transfer. At large distances the experimental slopes are generally in 
good agreement with the predictions of this model for one-neutron transfer. 
However, deviations from the expected ratio of two have been observed in 
two-nucleon transfer reactions, which in the literature are referred to as "slope 
anomalies" 0. Furthermore, some experiments show an energy dependence 
of the slope ||. 

In previous work [|], |5|, |J, we were able to explain available proton and 
neutron transfer data, including the slope anomaly, by considering the con- 
tribution of the two trajectories that lead to a given scattering angle (see 
Section II) . Thus far, we have applied this model exclusively to reactions for 
which structure effects in the transfer process are assumed unimportant. 

In the present work we investigate the transfer process for the case of 
vibrational nuclei, in particular one- and two- proton stripping reactions, at 
energies close to the barrier. In Section II we provide a brief review of the 
semiclassical theory of transfer in the case of structureless spherical nuclei, 
studied before, and extend it to the vibrational nuclei considered here. The 
results of the calculations performed with the model are presented in Section 
III. In the last section we draw the main conclusions and suggestions for 
additional work. 



P tr oc sin (0/2) e~ 2KDR 



(1) 



with 




(2) 
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II. THEORY 



a. Structureless nuclei 

In the framework of this semiclassical model 0, [5], the trajectories of 
the participant ions are determined taking into account both the Coulomb 
and the nuclear part of the nucleus-nucleus interaction. This brings, as a 
consequence, the possibility that more than one trajectory contributes to 
a given scattering angle of the outgoing particle. Due to absorption usually 
only two of them contributes to the transfer reactions. For the nuclear optical 
potential we adopt a Woods-Saxon shape with radius and strength calculated 
as in Ref. 0: 

R = R p + R t + 0.29 fm (3) 

with 

Ri = (1.2334 1/3 - 0.984~ 1/3 ) f m i = p,t (4) 

and 



V = ldTrjRa MeV 



(5) 



with 



7 = 0.95 



1 - U 



n p - z; 

An 



A f 



and 



R 



RpRt 
Rp + Rt 



MeV fm" 



(6) 



(7) 



where a is the diffuseness, and A i} N i} and Zj are the mass, neutron, and 
atomic numbers of the nucleus i (i = p for projectile, i — t for target), 
respectively. 

The probability amplitude for survival from absorption due to the imagi- 
nary part of the nucleus-nucleus optical potential, W(r), is calculated by the 
expression H 



a abs = exp 



W(t)dt 



(8) 
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in which W(t) = W{r{t)). 

We denote by U(r) the potential which acts over the transferred particle, 
Rb the position where this potential barrier reaches its maximum, U B = 
U(Rb), and B.E. the binding energy of the particle in the donor nucleus. 
The probability for tunneling is determined, when Ub + B.E. > 0, by the 
WKB approximation 

Ptun = Wtunl 2 = (l + e 5 )" 1 (9) 

in which 



R 2 



1/2 



5 = 2/ -£(U(r)+B.E.) dr. (10) 
JRi L n 

In the region Ub + B.E. < 0, the potential barrier can be approximated by 
an inverted parabola, allowing us the use of the analytic expression of Hill 
and Wheeler || 

P tun = [1 + exp (2ir/hu{U B + B.E.))]- 1 (11) 

with 

/ h 2 d 2 u(R B )\ 1/2 

hu=(-- -V^ • (12) 



/i dr z 



We have taken 



U(r) = U 1 (r) + U 2 (D-r), (13) 

U 1 (r) = U Cj (r) + U Nj (r), (14) 

where the subscripts 1 and 2 refer to the donor and acceptor cores, respec- 
tively, D is the distance of closest approach between them, and r is the spatial 
coordinate of the transferred particle with respect to the donor core. Uq is 
the Coulomb potential and U^. the nuclear potential generated by the core 
j over the particle. The Coulomb potential was taken as that generated 

1 /3 

by a charged sphere of radius 1.25 Aj fm acting over a cluster with charge 
Z c i and the nuclear part as a Saxon- Woods potential with radius parameter 
ro = 1.2 fm, diffuseness a u = 0.63 fm and depths 
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U Q] = Z d V 0] (+l) + N d V 0] (-l) 

where 



(15) 



51 - 33r, 



iV, - Z, 



3 ' 

•J 



MeV . (16) 
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is the nuclear potential generated by the core j with mass number Aj, charge 
Zj and neutron number Nj = Aj — Zj acting over the neutrons (r 2 = +1) 
and the protons (r 2 = -1) of the transferred cluster with mass number A d , 
charge Z d , and neutron number N d = A d — Z d . 

As detailed in Ref . |3J] , except in cases in which the measurements are done 
with high angular resolution, the transfer probability can be approximated 
by the incoherent sum of contributions by each trajectory leading to a given 
scattering angle 

P tr (e)^Y, P tun{0) \a abs (0) | 2 • (17) 

This expression is employed in the calculations presented in this work, 
b. Vibrational nuclei 

For a transfer process in which one of the participant ions is a vibrational 
nucleus (the target nucleus, to fix ideas) the above mentioned model can be 
extended in the following way. Assuming a quadrupole vibrational mode, we 
parametrize the radius of the target as 

Rt(a 20 ) = Rt(l + a 2Q Y 20 (9)). (18) 
The internal Hamiltonian of the target may be written as |L(J 

H int = -{B 2 \a 20 \ 2 + C 2 \a 20 \ 2 ). (19) 
from which the zero-point amplitude is 

<<wj (20) 

The target radius Rt is given by Eq. and E 2+ is the transition energy 
between the first excited and the ground state. 
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In this kind of reaction the vibration is very slow in comparison to the 
translational motion of the projectile and the tunneling process is dominated 
by those trajectories for which the distance between the surfaces of projectile 
and target at the point of closest approach are smallest. Therefore the most 
relevant axis of vibration is that directed from this point to the center of the 
target. In this spirit we substitute the spherical harmonic in Eq. ([18]) by its 



maximun value \Jhj 4-7T. 

With these prescriptions we now calculate the transfer probability as 
described for structureless nuclei, by substitution of Eqs. (§) and (g) or (|TT|) 
into Eq. ([[7]). The radius of the target must be replaced by i? 4 (a 20 ), which 
implies an additional degree of freedom in the calculation. Then, we solve 
the classical equations of motion using the Hamiltonian 

H = % + H int + V (r, a M ) (21) 

where T r is the kinetic energy of the relative motion and V(r, 020) — Vc + 
Vn includes the potential energy of the relative motion and the interaction 
between this and the vibrational mode. We approximate 

Vc = Z p Z t e> ( 1 _ + l R P R ^o) \ (22) 



for the Coulomb interaction and write 

VN = 1 + e [r-R P -Rt(a 20 )]/a ( 23 ) 

for the nuclear part. In the calculation of the absorption, the imaginary part 
of the nuclear optical potential acting between the incident ions is taken as 

WN = 1 + e [r~R p ~R t (a 20 )]/a ^ 

We integrate the resulting coupled differential equations with the initial 
condition a(t = — oo) = a^ cos0, where 0, the initial phase, can take any 
value in [0, 27r] with equal probability. For each phase the distance of 
closest approach and the value of a 2 o a t this time are calculated, in order to 
calculate the transfer probability for this trajectory. For a given scattering 
angle the contributions are summed as in Eq. flTT| ) and the final transfer 
probability is obtained by averaging over all phases. 
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III. RESULTS AND DISCUSSION 

We apply this model to the vibrational nucleus 98 Mo, studied by means 
of the 37 C1 + 98 Mo reaction. Data for reactions with Mo isotopes, including 
this system, at energies close to the barrier, were recently measured by the 



SUNY group Rill , [L2|. It was found that, in the case of two-proton stripping 
reactions, there is an abrupt change in slope at approximately the energy 
corresponding to the Coulomb barrier (Ej a & ~ 117 MeV). A calculation as 
described in Section II. a, i.e. without considering the vibrational character 
of the Mo nuclei, yields a value for the energy at which the slope changes 
higher than the observed one. 

In Fig. 1 we show the experimental transfer probabilities divided by 
sin(0 c . m ./2) for one-proton (circles) and two-proton (squares) stripping re- 
actions in the 37 CI + 98 Mo system ||12j|. Also shown are the theoretical 



results calculated as described in Section II. b and normalized to the data 
(full lines). We used the parameters of Temmer and Heydenburg [[Uj for 
98 Mo, E 2+ = 0.786 MeV, C 2 = 70 MeV. For the optical potential we take a 
diffuseness a = 0.7 fm, and a strength for the imaginary part Wq = 70 MeV. 
As shown in Fig. 1, a remarkably good agreement between the calculated 
and the experimental points is obtained at large values of D^ u th- in this 
approach D Ruth is only a parametrization of the scattering angle ||, ^ |6|]. 

We consider particularly interesting that our calculation reproduces the 
abrupt change in slope observed at bombarding energies around the Coulomb 
barrier. To understand the origin of this effect we should remember that, for 
a given deflection angle, there are two trajectories contributing to the transfer 
probability 0. One is essentially a Rutherford trajectory while the other feels 
more strongly the nuclear optical potential. At lower energies the dominant 
contribution to the transfer probability is the Coulomb trajectory, which 
gives the energy independent slope of Eq. (1), while at higher energies the 
nuclear trajectory dominates. As shown in Ref. || there is a rapid transition 
between these two regimes, which, for a structureless nucleus takes place 
at an energy well above the Coulomb barrier. In the case of a vibrational 
nucleus, the effective radius increases, which decreases the effective barrier, 
bringing, as a consequence, the transition in slope to lower energies. 

IV. CONCLUSIONS 

An extension of a simple semiclassical model including some nuclear struc- 
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ture effects, quadrupole vibrations to be specific, was proved to be quite suc- 
cesful. Applications to other degrees of freedom, such as octupole vibrations, 
can be implemented in essentially the same way. One could, in a similar 
spirit, include other nuclear properties, such as deformation. The practical 
problem with the nuclear deformation case is that the trajectories depend 
on the orientation of the deformation axis given by two Euler angles. The 
averaging process should then be performed over two parameters, the initial 
values of the two Euler angles and, consequently, one would need to consider 
a very large number of trajectories contributing to a given scattering angle, 
thus making the problem much more difficult to treat. 
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Figure Caption 

Fig. 1 One- and two- proton transfer probability, divided by sin(0 c . m ./2), as 
a function of D^ ut ^ for the 37 CI + 98 Mo reaction at three laboratory en- 
ergies. Symbols represent the experimental data of Ref. ||12|| , and lines 



are the theoretical calculations described in this work, and normalized 
to the data. 
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